Cooper pair dispersion relation for weak to strong coupling 
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Cooper pairing in two dimensions is analyzed with a set of renormalized equations to determine 
its binding energy for any fermion number density and all coupling assuming a generic pairwisc 
residual interfermion interaction. Also considered are Cooper pairs (CPs) with nonzero center-of- 
mass momentum (CMM) — usually neglected in BCS theory — and their binding energy is expanded 
analytically in powers of the CMM up to quadratic terms. A Fermi-sea-dependent linear term in the 
CMM dominates the pair excitation energy in weak coupling (also called the BCS regime) while the 
more familiar quadratic term prevails in strong coupling (the Bose regime). The crossover, though 
strictly unrelated to BCS theory per se, is studied numerically as it is expected to play a central 
role in a model of superconductivity as a Bose-Einstein condensation of CPs where the transition 
temperature vanishes for all dimensionality d < 2 for quadratic dispersion, but is nonzero for all 
d > 1 for linear dispersion. 
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The original Cooper pair (CP) problem |jj in two (2D) 
and three (3D) dimensions possesses ultraviolet diver- 
gences in momentum space that are usually removed 
via interactions regularized with large-momentum cut- 
offs [||. One such regularized potential is the BCS model 
interaction which is of great practical use in studying 
Cooper pairing Q and superconductivity |J. Although 
there are controversies over the precise pairing mecha- 
nism, and thus over the microscopic Hamiltonian appro- 
priate for high-T c superconductors, some of the proper- 
ties of these materials have been explained satisfactorily 
within a BCS-Bose crossover picture Q-Q] via a renor- 
malized BCS theory for a short-range interaction. In 
the weak-coupling limit of the BCS-Bose crossover de- 
scription one recovers the pure mean-field BCS theory 
of weakly-bound, severely-overlapping CPs. For strong 
coupling (and/or low density) well-separated, nonover- 
lapping (so-called "local") pairs appear [Q in what is 
known as the Bose regime. It is of interest to detail how 
renormalized Cooper pairing itself evolves independently 
of the BCS-Bose crossover picture in order to then discuss 
the possible Bose-Einstein (BE) condensation (BEC) of 
such pairs. We address this here in a single-CP pic- 
ture, while considering also the important case (gener- 
ally neglected in BCS theory) of non-zero center-of-mass- 
momentum (CMM) CPs that are expected to play a sig- 
nificant role in BE condensates at higher temperatures. 

In this Report we derive a renormalized Cooper equa- 
tion for a pair of fermions interacting via either a zero- 
or a finite-range interaction. We find an analytic expres- 
sion for the CP excitation energy up to terms quadratic 
in the CMM which is valid for any coupling. For weak 
coupling only the linear term dominates, as it also does 
for the BCS model interaction || . The linear term was 
mentioned for 3D as far back as 1964 (Ref. M p. 33). 



For strong coupling we now find that the quadratic term 
dominates and is just the kinetic energy of the strongly- 
bound composite pair moving in vacuum. 

The CP dispersion relation enters into each summand 
in the BE distribution function of the boson number 
equation from which the critical temperature T c of BEC 
of CPs is extracted. The linear dependence on the CMM 
of CP binding for weak coupling leads to novel transi- 
tion properties even in a heuristic BEC picture of su- 
perconductivity as BE-condensing CPs. It is well- 
known that BEC is possible only for dimensionalities 
d > 2 for usual nonrelativistic bosons with quadratic 
dispersion; this limitation reappears in virtually all BEC 
schemes thus far applied to describe superconductivity 
(Refs. [[ll 12 1 among others). But for bosons with a lin- 
ear dispersion relation as found here in weak coupling, 
BEC can now occur for all cZ > 1. Here we discuss 
the CP dispersion relation only in 2D. We have also per- 
formed a similar analysis in 3D and obtained the linear 
to quadratic crossover in the dispersion relation. 

Consider a two-fermion system in the CM frame, with 
each fermion of mass m , interacting via the purely at- 
tractive short-range separable potential pT[ 



V pq = -v g p g q , 



(0.1) 



where vq > is the interaction strength and the (? p 's 
are the dimensionless form factors g = (1 + p 2 /Po)^ 1 ^ 2 1 
where the parameter po is the inverse range of the poten- 
tial so that, e.g., po — > oo implies g p — 1 and corresponds 
to the contact or delt a potential V(r) = — vqS(t). The 
interaction model (3.1) mimics a wide variety of possible 
dynamical mechanisms in superconductors: mediated by 
phonons, or plasmons, or excitons, or magnons, etc., or 
even a purely electronic interaction. In the first instance 
mentioned, two terms of the form (1) can simulate a 



1 



coulombic interfermion repulsion surrounded by a longer- 
ranged electron-phonon attraction. The momentum- 
space Schrodinger eigenvalue equation for a two-particle 
bound state in vacuum with binding energy B 2 > for 
interaction (3.1) is M 
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h 2 k 2 /r 



(0.2) 



where k is the wavenumber in the CM frame and 
h 2 k 2 /2m the single-particle energy. 

On the other hand, the CP equation for two fermions 
above the Fermi surface with momenta wavevectors ki 
and k 2 (and arbitrary CMM wavevector K = ki + k 2 ) is 
given by 
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Cfe — — ^ VkqC q , 



(0.3) 



where k = |(k 



1 — k 2 ) is the relative momentum, Ek = 
2Ep — Ak the total pair energy, Ak > the CP bind- 
ing energy, C q = (q \ \&) its wave function in mo- 
mentum space, and the prime on the summation im- 
plies restriction to states above the Fermi surf ace: viz., 
|k ± K/2| > kp. For the separable interaction ( p~l| ) Eq. 
( |0.3|) becomes 
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(0.4) 



Although the summand in Eq. ((^4) is angle- 
independent, the restriction on the sum arising from the 
filled Fermi sea is a function of the relative wave vector k, 
and therefore angle- dependent. Th e po tentia l str ength vq 
can be eliminated between Eqs. ([l^) and ( |0.4[ ) leading 
to the renormalized CP equation 
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B 2 + H 2 k 2 /m 



9l 



h 2 k 2 /m + A K - 2E F + h 2 K 2 /Am 



(0.5) 



Instead of the arbitrary cutoff usu ally employed in deal- 
ing with delta interactions, in Eq. ( |0.5| ) we rely on physi- 
cal "observables" for the sake of renormalization, viz., the 
ground- stat e binding energy B 2 in vacuum. The sums 
in Eq. (3J3) can be transformed to integrals; the restric- 
tion in the second term arising from the filled Fermi sea 
leads to two different expressions depending on whether 
K = K/kp is < 2 or > 2, as discussed in the Ap- 
pendix. Letting all variables be dimensionless by ex- 
pressing them either in units of the Fermi wavenum- 
ber kp or of the Fermi energy Ep = h 2 k 2 7 /2m 1 viz., 
£ = k/kp, B 2 = B 2 /E F , A K = A K /E F , etc., we define 
o? K = l-A K /2- K 2 /4 ee -0k, and 9 the angle between 

wavevectors k and K so that £o($) = y 1 — K 2 sin 2 6/4 + 



Kcos9/2 and &(0) ee -yT - K 2 sin 2 9/4 + Kcos9/2. 
For a zero-range interaction, = 1, after some algebra 
one gets 
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d9He {9)-a 2 K ] = \\u^. 



K < 2, 
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(0.6) 

K > 2, 
(0.7) 



where 9 = arcsin(2/AT) < 7r/2. For K = Q only Eq. 



( p.6[ ) applies, in which case ^o(^) = 1, a z K = = 1 — 
Ao/2 and we obtain the surprising result Aq = B 2 , i.e., 
for an attractive delta interaction the vacuum and CP 
binding energies for zero CMM coincide for all coupling, a 
result apparently first obtained in Ref . |^] . For a nonzero 
CMM the CP binding energies can be calculated 
from Eqs. (0.6) and ( 0.7). For K ^ a minimum 



threshold value of B 2 /Ep is found to be required to bind 
a CP. 




Equations (0.6) and ( p.7[ ) can be solved numerically 
for the CP binding A K for any CMM. For small CMM 
only Eq. ( |0.6| ) is relevant; and this equation for small but 
non-zero K and for K — can be subtracted one from 
the other. This gives the small-CMM expansion valid for 
any coupling B 2 /Ep > 0, 



0(K 3 ), (0.8) 



where a nonnegative CP excitation energy sk has been 
defined, and the Fermi veloci ty vp comes from Ep /kp = 
hvp/2. The leading term in jo^) is linear in the CMM, 
followed by a quadratic term. The linear CP dispersion 
term should not be confused with that of the many-body 
(collective) excitation spectrum in weak coupling. Only 
CPs can undergo BEC while bosonic "excitations" (or 
modes or phonons) cannot since the former are fixed in 
number while the latter are not. Indeed, the particle- 
hole [sometimes called the Anderson-Bogoliubov-Higgs 
(ABH)] modes of excitation energy hvpK/Vd [[l3| in d 
dimensions in the zero coupling limit are larger than 
the weak-coupling CP dispersion energies (2/n)hvpK 
and \hv F K (Ref. § p. 33) in 2D and 3D @, respec- 
tively, while in ID they happen to coincide — in spite of 
the fact that CPs and ABH-like modes are physically 
dist inct entities. The coefficient of the quadratic term 
in (|J) changes sign at B 2 /E F = A /E F ~ 0.379, 
as one goes from weak (B 2 = Ao <§; Ep) to strong 
(B 2 = A ^> Ep) coupling. If vp (or Ep) — > explic- 
itly (dilute limit) the first two terms of Eq. ( p.8[ ) reduce 
simply to 
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(0.9) 



for any coupling. This is clearly just the familiar non- 
relativistic kinetic energy in vacuum of the composite 
(so-called "local") pair of mass 2m and CMM K. The 
same result ( fL9| ) is also found to hold in 3D. 

Figure 1 shows exact numerical results for the zero- 
range potential (gk = 1) for different couplings of a CP 
excitation energy ek/A as function of CMM K/kp, 
both dimensionless. We note that the CPs break up when- 
ever Ajf tu rns from positive to negative, i.e., vanishes, or 
by Eq. (3.£) when e k/A q — 1. These points are marked 
in the figure by dots. In add ition to t he exact results 
obtained by solving Eqs. ( |0.6|) and ( |0.7|) , we also exhibit 
the results for the linear app roximation [first term on the 
right-hand side of Eq. (3.8), dot-dashed lines, virtually 
coinciding with the exact curve for all B^/Ep < 0.1] as 
well as for the quadrati c ap proximation (dashed parabo- 
las) as given by Eq. (D.9) for stronger couplings. For 
weak enough coupling or large enough fermion density 
at any nonzero coupling the exact dispersion relation is 
virtually linear — in spi te o f the divergence of the isolated 
quadratic term in Eq. ( pTq ) as B^/Ep — > 0. As c oupling 
is increased the quadratic dispersion relation (0.9) slowly 
begins to dominate. The crossover from a linear to a 
quadratic dispersion relation manifests itself by a change 
in curvature from concave down to concave up — these 
two regions being separated by an inflection point that 
moves down towards the origin as coupling is increased 
to infinity. 

Figure 2 shows the CP excitation energy ejf/A as 
a function of CMM K/kp calculated for the finite-range 
interaction form factor g p — (l+p 2 /p§) -1 / 2 withpo = kp 
for weak to moderate coupling. To compare, we also plot 
the zero-range result as well as the linear rela tion given 
by the first term on the right-hand side of Eq. ( |0.8| ). The 
finite-range curves are closer to the corresponding zero- 
range ones if labeled by Aq/ Ep instead of by B2/EP , as 
was done with all four sets of curves. 

Figure 3 exhibits Sk/A q as a function of Kjkp for 
the finite-range interaction with po = kp for stronger 
couplings. In this case there is no special advantage in 
labeling the dispersion curves by Ao so Bi was used with 
results for B2/EP = 3, 10 and 20 shown. In the zero- 
range case the curves gradually tend to the quadratic 
form as B2 increases. For finite-range, po — kp, the 
curves develop a maximum followed by a minimum with 
a point of inflection in between. The slope at the point of 
inflection is now negative. Although each curve tends to 
a quadratic form for large enough K/kp, they are quite 
different from it for small K/kp. These "looped" dis- 
persion curves are reminiscent of the "roton" excitation 
spectrum |l4|] in liquid 4 He. 

To summarize, the single CP problem with non-zero 
CMM is tracked as it evolves in varying the interfermion 
short-range pair interaction from weak to strong or vary- 
ing fermion density from high to low, respectively, for 



any fixed nonzero coupling. The CP excitation energy 
is exhibited as a function of its CMM. For weak cou- 
pling, the excitation energy is a linear dispersion relation 
in the CMM, and changes very gradually to a quadratic 
relation as coupling increases. For a zero-range pair in- 
teraction in the dispersion curve one typically has a point 
of inflection with a positive slope separating a region of 
concave-down curvature for small CMM from a region 
of concave-up curvature for large CMM. For finite-range 
pair interactions of sufficiently long range the slope at 
the point of inflection changes from positive to zero and 
eventually becomes negative. This leads to maxima and 
"roton-like" minima in the CP dispersion curves. These 
results will play a critical role in a model of supercon- 
ductivity based on BE condensation of CPs as they will 
yield, even in 2D as in the cuprates , BEC transition tem- 
peratures T c that interpolate between nonzero values in 
weak coupling with a linear CP dispersion relation down 
to the expected T c = value in strong coupling with a 
quadratic relation. 

Appendix: The restriction that both particles lie 



above the Fermi sea in Eq. (0.5) can be written as 



{k/kp ±K/2k F f - 1 = £ 2 ±£Kcos9 + K 2 /4- 1 > 0, 

(0.10) 

where £ = k/kp and K = K/kp. The equality leads to 
two pairs of roots in £, say £1,2 = — a±b and £3,4 — a±b, 



where a = (K/2)cos9, b = yi- (K 2 /4) sin 2 9, and 9 
the angle between k and K. 

For K < 2, b > a, one root of the two pairs is pos- 
itive and the other negative. Thus, Eq. ( |0.10| ) can be 
satisfied provided that £ > £1,^2, £3, £4, or specifically, if 
£>£ (#) = a + fc. 

For K > 2 and 9 > 0q = arcsin(2/ii"), b becomes 
imaginary and Eq. ( 0.1 0| ) is satisfied for all £. Therefore, 
there is no restriction in the integration over £. However, 
for K > 2 and 9 < 9q, b < a the pair of roots £1,2 are both 
negative while the pair £3^ are both positive (with £3 > 
£4). Consequently, in both cases Eq. ( p. 10 ) is satisfied 
only if £ is in the interval [0, £q(0) = a — b], and in the 
interval [£o(#),oo] , respective ly. Using these restrictions 
on the £ inte gratio n in Eq. ( |0.5[ ) one eventually arrives 
at Eqs. (|0j) and ( |b~7|) . 
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Fig. 1. Dimensionless CP excitation energy e^/Ao 



vs K/kp , calculated from Eqs. (0.6) and (3.7) for differ- 
ent couplings B2/EF, full curves. The dot-dashed line 
is the linear approximation (virtually coincident with the 
exact curve for B^jEp < 0. 1) w hile the dashed curve is 



the quadratic term of Eq. (0.9). Dots denote values of 



CMM wavenumber where the CP breaks up, i.e., where 
A/f = 0. 

Fig. 2. Same as Fig. 1 but for couplings expressed 
as Ao/Ep. The dot-dashed line is the linear approxima- 
tion; the dashed curve is the result for the finite-range 
interaction pq — hp, and the full curve is the zero-range 
result. For the finite-range potential Ao/ Ep = 0.1, 0.5, 
1.0 and 2.0 correspond to B 2 /E F = 0.469, 1.4, 2.45 and 
4, respectively. Dots and squares mark values of CMM 
wavenumber where the CP breaks up. 

Fig. 3. Same as Fig. 1 but including also finite range 
at stronger couplings B2/EP. The full curve is the exact 
zero-range result; the short-dashed one the quadratic ap- 
proximation; the long-dashed one the exact finite-range 
result with p = hp. Each set of three curves is labeled 
by different values of B 2 /Ep. 
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